Assessment of urban ecosystem health and its influencing factors: a case study of Zibo City, China

Urban ecosystem health is the foundation of sustainable urban development. It is important to know the health status of urban ecosystem and its influencing factors for formulating scientific urban development planning. Taking Zibo city as the study area, the indicators were selected from five aspects: ecosystem vigor, structure, resilience, service function and population health to establish an assessment index system of urban ecosystem health. The health level of urban ecosystem was assessed, and its changing trend was analyzed from 2006 to 2018 in Zibo. Furthermore, obstacle degree analysis and sensitivity analysis were used to quantitatively analyze the main obstacle factors and sensitivity factors affecting urban ecosystem health, so as to provide references for improving urban ecosystem health. The results showed that the health level of urban ecosystem in Zibo showed an upward trend from 2006 to 2018. The poor structure and ecological environment quality were the main obstacle factors to urban ecosystem health. The impact of changes in a single indicator on urban ecosystem health gradually decreased, but the sensitivity index of indicators had obvious differences. Urban ecosystem health was sensitive to changes in ecosystem structure and resilience. In the future, Zibo should strengthen ecological construction, optimize the industrial structure, and develop green economy to promote urban ecosystem healthy.


Study area and data sources
Zibo city lies between 117° 32′-118° 31′ E and 35° 55′-37° 17′ N, with a total area of 5965 km 2 .In terms of topography, mountains are distributed in the east, south and west of the region, and the central part is lower and northward-inclined.The terrain is high in the south and low in the north, and the height difference between north and south is more than 1000 m.Centrally located in Shandong province, Zibo is not only an important transportation hub in the province with developed transportation, but also an important petrochemical, pharmaceutical production base and building materials production area in China with rich mineral resources.Due to its advantages of location and resource base, Zibo has become an important industrial city.However, the regional ecosystem health is facing the adverse effects of traditional industrial production models and industrial structures.Therefore, the urgent task for researchers and policymakers is to accelerate the transformation of urban economies and promote urban ecosystem health.
Taking Zibo city, Shandong province as the study area, the indicator data were collected from the municipal districts of Zibo.These data were sourced from China City Statistical Yearbook, Shandong Statistical Yearbook, Zibo Statistical Yearbook and Zibo Statistical Bulletin on National Economic and Social Development from 2007 to 2019.Among them, 7 indicators were obtained from China City Statistical Yearbook, including Per capita GDP, GDP growth rate, proportion of tertiary industry in GDP, urban sewage treatment rate, harmless treatment rate of domestic waste, comprehensive utilization rate of industrial solid waste and urban road area per capita, energy consumption per GDP and the proportion of environmental protection investment in GDP were retrieved from Shandong Statistical Yearbook, days with air quality of grade II and above was extracted from Zibo Statistical Bulletin on National Economic and Social Development, and the remaining 10 indictors were sourced from Zibo Statistical Yearbook.

Establishing the assessment index system
According to the meaning of urban ecosystem health, the vigor-structure-resilience-function-population health framework was adopted to construct the assessment index system in this study, which has been well applied in relevant research 26,31 .From previous research 11,22,23,45,46 , it can be seen that due to the differences in natural conditions, economic development level, industrial structure and development goals of different cities, although the selected evaluation indicators are different, this evaluation framework can reflect the productivity, structural complexity, stability, ability to provide services and maintain human health in urban ecosystem.We first sorted out and analyzed the evaluation indicators in the existing literature, and preliminarily selected some evaluation indicators.Then, according to the development characteristics of Zibo city, the preliminarily selected indicators were further analyzed and screened.
Life expectancy and under 5 mortality rate were selected as indicators of population health in previous studies.Since under 5 mortality rate data were incomplete in this study period and it only reflected the health status of children, and it is difficult to obtain data on life expectancy, we selected the mortality rate to reflect the overall health status of urban residents.From the existing indicator system, it can be found that there was a lack of indicators that reflect the medical health level of urban residents.In this study, health care expenditure of urban residents was added to reflect the ability of urban ecosystem to maintain human health.In addition, considering the characteristics of industrial city and the actual situation of Zibo city, some indicators were appropriately added.Days with air quality of grade II and above was selected to reflect air quality, proportion of high-tech industries in GDP was selected to reflect the industrial structure, two indicators, harmless treatment rate of domestic waste and proportion of environmental protection investment in GDP were selected to reflect the self-regulation ability of urban ecosystem under human management activities.Finally, evaluation indicators were selected from five aspects: ecosystem vigor, structure, resilience, service function and population health to construct the assessment index system of urban ecosystem health in Zibo (Table 1).(a) Vigor, which reflects a city's economic vitality and green production efficiency.We selected three indicators: Per capital GDP (x 1 ), GDP growth rate (x 2 ) and Energy consumptions per GDP (x 3 ).(b) Structure, which reveals the natural, economic and social structure in urban areas.We selected five indicators: Proportion of tertiary industry in GDP (x 4 ), Proportion of high-tech industries in GDP (x 5 ), Registered urban unemployment rate (x 6 ), Urban population density (x 7 ) and Green coverage rate in built-up area (x 8 ).(c) Resilience, which reflects the self-regulating capacity of urban ecosystem.Since urban ecosystem is a highly artificial ecosystem, its self-regulating capacity is closely related to human management activities.Thus, Urban sewage treatment rate (x 9 ), Comprehensive utilization rate of industrial solid waste (x 10 ), Harmless treatment rate of domestic waste (x 11 ), and Proportion of environmental protection investment in GDP (x 12 ) were selected.(d) Service function, which reflects the ability of urban ecosystem to provide production and living services for urban residents.Per capita housing construction area (x 13 ), Urban road area per capita (x 14 ), Number of hospital beds per 10,000 people (x 15 ), Days with air quality of grade II and above (x 16 ), and Park green space per capita (x 17 ) were selected.(e) Population health, which reflects the people-oriented characteristics of urban ecosystem.Engel's coefficient of urban residents (x 18 ), mortality rate (x 19 ), and Health care expenditure of urban residents (x 20 ) were selected to reflect the living standard and health status of urban residents.

Data standardization
In order to eliminate the influence of different units, features, and orders of magnitude of indicators, the original data should be firstly standardized.According to the relationship between the indicator and the evaluation goal, the indicator can be divided into three types, positive indicators, negative indicators and appropriate indicators.(2011-2015) and "13th Five-Year" Plan (2016-2020) of Zibo city, 10% is taken as the ideal GDP growth rate.The scale effect brought by urban population agglomeration is beneficial to urban development, but excessive population density will also lead to negative effects such as "urban disease".Therefore, urban population density should be controlled reasonably.In China, urban land area is allocated according to urban population in urban development planning 47 .According to the standard of 100 m 2 per capita, the urban population density should reach 10,000 people per square kilometer.In addition, considering the scarcity of land resources and environmental protection, the ideal values of appropriate indicators are determined by referring to the recommended values of ecological city, healthy city, garden city and environmental protection city at home and abroad to improve the utilization efficiency of public resources and rationally adjust land use structure (Table 1).
The values of various indicators are standardized to make them between 0 and 1.The formula is as follows: where x ij is the original value of indictor j in year i; x ' ij is the normalized value of indictor j in year i; x jmax is the maximum value of indictor j in the study period; x jmin is the minimum value of indictor j in the study period; and x * j is the ideal value of indictor j.

Weight calculation
The entropy weight method is used to calculate the weight of each indictor, and the formula is as follows: Calculating the proportion of indictor j, R ij

Calculating the entropy of indictor j, e j
Calculating the difference coefficient of indictor j, g j Calculating the weight of indictor j, w j where x ' ij is the normalized value of indictor j in year i; m is the number of the assessed years; and n is the number of indicators.
Since the entropy weight method uses the concepts of entropy and logarithm, subject to relevant rules, negative value cannot be involved in the calculation.Then, when R ij = 0, 0.00001 is used instead.The calculation results of indictor weight are shown in Table 1.

Urban ecosystem health assessment model
The multi-factor comprehensive evaluation model is adopted to calculate the urban ecosystem health index.The formula is as follows: (1) where F is the urban ecosystem health index; x ' ij is the normalized value of indictor j; w j is the weight of indictor j; w i is the weight of assessment factor i; m is the number of assessment factors; and n is the number of indicators in the assessment factor i.

Rank correlation analysis method
Spearman rank correlation analysis is used to analyze the changing trend of urban ecosystem health index.The calculation can be expressed as: where r s is the Spearman rank correlation coefficient; x i and y i are the serial number of urban ecosystem health index and year from small to large, respectively, namely rank order; d i is the difference between x i and y i ; and N is the number of samples.

Obstacle degree method
By introducing contribution index, deviation degree index, and obstacle degree index, the obstacle degree model is used to diagnose the obstacle factors that affect urban ecosystem health.The formula is as follows: where V j is the contribution degree of indictor j; P j is the deviation degree of indictor j; A j is the obstacle degree index of indictor j; w j is the weight of indictor j; w i is the weight of factor i to which indictor j belongs; x j is the value of indictor j; and n is the number of indictors.

Sensitivity analysis method
Referring to relevant research [48][49][50] , the elasticity analysis model is used to calculate the impact intensity of indictor changes on urban ecosystem health.The sensitivity index is defined as the ratio of the change rate of urban ecosystem health index to that of the indictor, that is, the change degree of urban ecosystem health index caused by 1% of the indictor change.The larger the sensitivity index, the more sensitive the urban ecosystem health index is to the change of indictor value, and the greater the impact of indicator change on urban ecosystem health.Conversely, the less sensitive it is, the smaller the impact of indicator change on urban ecosystem health.The formula is as follows: where CS is the sensitivity index; F m and F n are the initial index and the adjusted index of urban ecosystem health, respectively; x m and x n are the initial value and the adjusted value of indicator, respectively.

Analysis of urban ecosystem health change in Zibo
Overall change characteristics of urban ecosystem health On a whole, the health level of urban ecosystem in Zibo showed a fluctuating upward trend (Table 2).The health index of urban ecosystem increased from 0.329 in 2006 to 0.897 in 2018, which increased by 2.73 times.( 9) www.nature.com/scientificreports/At different stages of socio-economic development, obvious differences emerged in the health status and its changes of urban ecosystem in Zibo.From 2006 to 2010, the urban ecosystem health index was low, with an average value of 0.383, but showed a rapid improvement, with an average annual increase of 13.2%, which was consistent with the development characteristics of Zibo.On the one hand, as an industrial city, Zibo had a high proportion of the secondary industry, but a relatively low percentage of high-tech industry, accounting for about 20% of the total output value of industries above designated size, resulting in higher energy consumption and waste output.Despite the rapid growth of GDP in this period, the comprehensive utilization rate of industrial solid waste and environmental protection investment were relatively low because of the weak technical level and weak environmental awareness, leading to the lower health level of urban ecosystem.On the other hand, guided by the policy of "energy conservation and emission reduction" put forward in the "11th Five-Year Plan" (2006-2010), the new industrial production concept was applied to gradually reduce energy consumption and improve the comprehensive utilization rate of waste.Simultaneously, with the gradual increase of environmental protection investment, the ecological environment improved significantly, and the health level of urban ecosystem continued to improve.The urban ecosystem index increased by 52.76% from 2006 to 2010.
From 2011 to 2018, the health level of urban ecosystem kept improving, but the rate of improvement tended to be gentle.From 2011 to 2015, the average urban ecosystem health index was 0.613, with an average annual increase of 8.13%.During this period, under the guidance of the policies of "transformation of Old-New-Driving-Forces" and "energy conservation and emission reduction", Zibo paid attention to urban ecological construction.Urban living environment was constantly improved by adjusting the industrial structure and increasing investment in environmental protection.Although the GDP growth showed a slowdown during the "12th Five-Year Plan" (2011-2015), the overall health level of urban ecosystem improved.From 2016 to 2018, the average urban ecosystem health index reached 0.828, with an annual increase of 7.39%, indicating that under the background of ecological civilization construction, Zibo no longer seeks economic development at the expense of ecological environment.In the early stage of "13th Five-Year Plan" (2016-2018), the health status of urban ecosystem reached a high level, and urban ecosystem was in a stable development state.

Change of subsystem health
From the perspective of five subsystems, vigor index, structure index, resilience index and population health index all changed greatly, while service function index fluctuated, but kept at a high level, indicating that Zibo paid attention to people's livelihood in urban construction and offered relatively perfect infrastructure and public service facilities.
From 2006 to 2010, vigor index and resilience index were both at a low level.Although economy developed rapidly during this period, due to the influence of traditional industrial production concepts and technologies, high energy consumption, insufficient investment in environmental protection and low comprehensive utilization rate of waste led to low vitality and resilience of urban ecosystem in Zibo.This reflected that the ecological environment problems caused by rapid economic development in this period were more prominent, and the development mode at the expense of the environment was not conducive to ecosystem healthy.After 2010, guided by the "energy conservation and emission reduction" policy and the concept of new industrial production, Zibo strengthened urban ecological construction.Through technological innovation and increased environmental protection input, energy consumption per GDP reduced by 6% annually, harmless treatment rate of domestic waste reached 100%, and urban sewage treatment rate and the comprehensive utilization rate of industrial solid waste increased by 10.46% and 7.37% respectively, thus enhancing the vitality and resilience of urban ecosystem.
The structure index showed a fluctuating upward trend.From 2006 to 2015, the structure index was low, which was manifested as the unreasonable industrial structure, low proportion of high-tech industry and the underdeveloped tertiary industry in Zibo.After 2015, with the implementation of the Old-New-Driving-Forces conversion policy, Zibo actively promoted the transformation and upgrading of old industries.By transforming old industrial enterprises, introducing high-tech industries and developing the tertiary industry, the industrial structure was constantly improved and the structure of urban ecosystem was also constantly optimized.
The population health index fluctuated greatly.Especially in 2012-2014, there was a significant decline, which may be due to the influence of policies.Since 2012, the GDP growth rate in Zibo began to decline, which affected residents' income and expenditure.In 2012, the Engel's coefficient of urban residents increased by 4% compared with 2011, and the health care expenditure of urban residents also decreased correspondingly, by 14.54% compared with 2011.From 2015 to 2018, as the economy recovered and the medical security system continued to improve, the population health index continued to improve, and people's living conditions showed a positive trend.

The change trend of urban ecosystem health
When N = 13, critical value W P = 0.56 (P < 0.05), and results of Spearman rank correlation analysis are showed in Table 3.The rank correlation coefficient of urban ecosystem health index is positive and higher than w p , www.nature.com/scientificreports/indicating that the overall health level of Zibo's urban ecosystem showed a significant improvement trend from 2006 to 2018.The variation trends of subsystem health are slightly different.The rank correlation coefficients of vigor index, structure index, resilience index and population health index are all positive and greater than w p , while the rank correlation coefficient of service function index is positive but less than w p .These results showed that the vigor, structure, resilience and population health of urban ecosystem in Zibo changed significantly and showed a significant improvement trend from 2006 to 2018, while the service function improved, but the improvement trend was not significant.

Obstacle factor analysis
On the basis of calculating the factor obstacle degree for each year, the average value of factor obstacle degree of each stage was taken according to the stage characteristics of urban ecosystem health in Zibo.The factor obstacle degree was arranged from large to small, and the factors with a cumulative obstacle degree of more than 80% were regarded as the main obstacle factors affecting urban ecosystem health.
The main obstacle factors affecting urban ecosystem health in Zibo were different in different periods, and the number of factors showed a decreasing trend (Table 4).Compared with 2006-2010, the main obstacle factors decreased from 10 to 9 in 2011 to 2015, among which, the obstacle degree of urban sewage treatment rate (x 9 ) and per capita GDP (x 1 ) decreased significantly and were no longer the main obstacle factors, while the obstacle degree of days with air quality of grade II and above (x 16 ) increased to be a new main obstacle factor, and the other eight main obstacle factors did not change.In general, during the period of "11th Five-Year Plan" (2006-2010) and "12th Five-year Plan" (2011-2015), the tertiary industry and high-tech industry accounted for a relatively low proportion in Zibo due to the influence of the industrial structure.In the process of rapid economic growth, high energy consumption in industrial production, insufficient investment in environmental protection and poor ecological environment were the constraints to urban ecosystem health.At the same time, the lack of expenditure on health care was also the main factor affecting population health in this period.
From 2016 to 2018, the main obstacle factors were further reduced to seven, and the obstacle factors also changed significantly.With the implementation of policies to conserve energy, reduce emissions, and accelerate transformation of Old-New-Driving-Forces, in the early stage of "13th Five-Year Plan" (2016-2018), the industrial structure of Zibo was constantly optimized, economy developed steadily, and environmental protection investment was gradually increased.The obstacle degree of four factors: energy consumption per GDP (x 3 ), the proportion of tertiary industry in GDP (x 4 ), the proportion of environmental protection investment in GDP (x 12 ) and Engel's coefficient of urban residents (x 18 ) were significantly reduced, which were no longer main obstacle factors.However, in the process of socio-economic transformation, the decline in economic growth affected the vitality of urban ecosystem, and GDP growth rate (x 2 ) became the main obstacle factor.In addition, air quality, proportion of high-tech industry, comprehensive utilization of industrial solid waste, employment and health care of urban residents, and urban population density were still the main obstacle factors affecting urban ecosystem health.
It is worth noting that the proportion of high-tech industry in GDP (x 5 ), registered urban unemployment rate (x 6 ), urban population density (x 7 ) and health care expenditure of urban residents (x 20 ) are the main obstacle factors in each period.The proportion of high-tech industry in GDP (x 5 ), registered urban unemployment rate (x 6 ), and urban population density (x 7 ) are all indicators reflecting the structure of urban ecosystem, indicating that poor structure is the main reason restricting urban ecosystem health in Zibo.Therefore, in the future, Zibo should continue to optimize the industrial structure, attract talents, increase employment rate, and improve residents' medical security.

Sensitivity analysis
In 2006, 2010, 2015 and 2018, the values of each indicator were adjusted up (or down) by 1% respectively while keeping the values of other indicators unchanged to analyze the change of urban ecosystem health index caused by 1% change of indictor value (Table 5).The results showed that only the sensitivity index of the proportion of www.nature.com/scientificreports/environmental protection investment in GDP (x 12 ) increased from 2006 to 2018, whereas the sensitivity indexes of other indicators showed a decreasing trend, indicating that increasing environmental protection input had an important impact on improving urban ecosystem health in Zibo.Meanwhile, it also showed that before 2010, when the indicators were at a low level, improving the quality of a single indictor could effectively improve the health status of urban ecosystem.Conversely, a decline in the quality of an indicator would also lead to a decline in urban ecosystem health.After 2010, with the continuous improvement of the structure and function of urban ecosystem, the impact of changes in a single indictor on ecosystem health was getting weaker and weaker, reflecting the improvement of self-regulation and anti-interference ability of urban ecosystem, and the continuous improvement of urban ecosystem health in Zibo.
From the results of sensitivity analysis, it can be seen that although the sensitivity of urban ecosystem health to changes of various indicators is decreasing in Zibo, there are still obvious differences.Based on the sensitivity index, the indicators can be divided into two categories, namely sensitivity indictor and non-sensitivity indictor.Since 2015, the sensitivity indexes of the proportion of tertiary industry in GDP (x 4 ), the proportion of high-tech industry in GDP (x 5 ), the registered urban unemployment rate (x 6 ), urban sewage treatment rate (x 9 ), comprehensive utilization rate of industrial solid waste (x 10 ), harmless treatment rate of domestic waste (x 11 ) and the Engel's coefficient of urban residents (x 18 ) are all greater than 0.2, obviously higher than other indicators.Changes of these indicators have a relatively large impact on urban ecosystem health, hence they are designated as sensitivity indicators.Other indicators with lower sensitivity index have a weak impact on urban ecosystem health, and are designated as non-sensitive indicators.Among the sensitivity indicators, except the Engel's coefficient of urban residents (x 18 ) which reflects the living standard of urban residents, the other six indicators reflect the structure and resilience of urban ecosystem, indicating that the structure adjustment and resilience change of ecosystem have a prominent impact on urban ecosystem health.According to the characteristics of industrial city, Zibo should strengthen urban ecological construction, improve the comprehensive waste treatment capacity, and at the same time continue to optimize the industrial structure, introduce high-tech industries, increase the employment opportunities for residents to effectively promote urban ecosystem healthy.

Discussion
From the perspectives of ecology and system theory, an important feature of urban ecosystem health is rational structure and perfect function 11 .In addition, human beings are an important component of urban ecosystem, and urban ecosystem should also provide services for human beings and maintain human health.In this study, the vigor-structure-resilience-function-population health assessment framework was constructed to comprehensively reflect the overall characteristics of the urban ecosystem from different dimensions.Referring to previous studies and based on the assessment framework, we screened and supplemented the indicators to construct a relatively complete assessment index system in this study.Some indictors were added, such as air quality indicator reflecting the characteristics of industrial cities, proportion of high-tech industries indicator reflecting the industrial structure, harmless treatment of domestic waste indicator and environmental protection investment indicator reflecting the self-regulation ability of urban ecosystem, and medical security indicator reflecting the medical health level of urban residents, which was the improvement of the index system.Based on the complexity of urban ecosystem, although the indicators have been screened and supplemented in this study, the index system still needs to be improved due to the limitations of data acquisition and quantification.In future research, more data sources and indicator types need to be considered to further improve the index system and improve the accuracy of assessment.Although the index system established in this study is not the most ideal, it can still provide references for city managers and decision-makers to maintain urban ecosystem health in Zibo.This study analyzed the influencing factors of urban ecosystem health from the aspects of obstacle degree and sensitivity, which can not only quantitatively reflect the influence degree of different factors, but also reflect the change trend of the influence degree of factors.This is conducive to identifying the key factors affecting the healthy development of urban ecosystems.To a certain extent, these research methods can make up for the shortcomings of existing research methods.
From 2006 to 2018, the health level of urban ecosystem continued to improve in Zibo, which was closely related to the policies of "energy conservation and emission reduction" and "transformation of Old-New-Driving-Forces". www.nature.com/scientificreports/As an industrial city, ecosystem structure and resilience are currently the main factors affecting urban ecosystem health in Zibo.In the future, policies will continue to play an important role in improving urban ecosystem health.Therefore, urban green development policies should be formulated to balance the relationship between economic development and environmental protection, so as to promote urban ecosystem health.Zibo should adopt environmental protection policies to strengthen urban ecological construction and enhance the resilience of urban ecosystem, such as formulating strict pollutant emission standards, formulating incentive policies for technological innovation, clean, regenerated energy utilization policies, and increasing environmental protection input.In the second place, Zibo should formulate industrial development policies and seek for a new economic development model.Developing the tertiary industry and increasing the proportion of high-tech industries will help Zibo optimize its industrial structure and promote regional economic development.In addition, the lower urban population density is an obstacle factor affecting urban ecosystem health in Zibo.Talents are important for urban development.Studies have shown that improving the quality of the population is conducive to promoting healthy urban development 21,26 .Zibo is not attractive enough for talents now.In order to promote urban development, Zibo should learn from the experiences of other cities and formulate talent introduction policies to attract all kinds of talents to settle down.Due to the impact of COVID-19 broke out in 2019, some indicator data changed greatly.In order to avoid the impact of data anomalies on the research results, the study period was selected from 2006 to 2018.The epidemic will have a negative impact on socio-economic development, which in turn will affect urban ecosystem health.The change characteristics of urban ecosystem health under the impact of the epidemic will be further studied in future research.

Conclusion
Assessment of urban ecosystem health can help to identify the health status and limiting factors in an urban ecosystem, which can provide references for governments to formulate urban management strategies.In this study, an index system was constructed based on the vigor-structure-resilience-service function-population health framework to assess the health level of urban ecosystem from 2006 to2018 in Zibo.Obstacle degree analysis and sensitivity analysis were used to quantitatively analyze the main obstacle factors and sensitivity factors affecting urban ecosystem health.The conclusions are as follows: Under the guidance of the new industrialization concept and the policy of transformation of Old-New-Driving-Forces, the overall health level of urban ecosystem showed a rising trend from 2006 to 2018 by optimizing the industrial structure and strengthening urban ecological construction, especially ecosystem vigor, structure, resilience and population health were significantly improved.Although the service function was not significantly improved, it continued to maintain a high level.
From 2006 to 2018, the number of main obstacle factors affecting urban ecosystem health showed a decreasing trend.The poor structure and ecological environment quality were still the main constraints on urban ecosystem health.
The results of sensitivity analysis showed that from 2006 to 2018, except for environmental protection investment, the sensitivity index of other indicators showed a decreasing trend, but the sensitivity of indicators was still significantly different.Ecosystem health was more sensitive to changes in structure indicators and resilience indicators.The impact of changes in a single indicator on urban ecosystem health gradually decreased, which was the result of gradual improvement of urban ecosystem, and also indicated that focusing only on one aspect of urban development could not effectively improve level of urban ecosystem health in the future.Therefore, promoting the coordination of socio-economic development and ecological environmental protection is a realistic choice to ensure urban ecosystem health in Zibo.

Table 1 .
Assessment index system of urban ecosystem health in Zibo.

factor Factor weight Assessment indictor Indictor weight Indictor feature (ideal value)
Positive indicator suggests that a higher value indicates a healthier urban ecosystem, while negative indicator suggests the opposite.Appropriate indicator means that only when the indicator reaches a specific optimal level does it contribute to urban ecosystem health.Among the appropriate indicators, GDP growth rate reflects the speed of urban economic development.From the perspective of sustainable development, rapid GDP growth is not always optimal, and moderate growth is more conducive to the healthy development of cities.During the period of "11th Five-Year Plan" (2006-2010), GDP grew by 14.21% annually in Zibo, higher than the national average.Since 2011, the strategy of sustainable development has been carried out, which no longer pursuing rapid growth in the economic quantity, but paying more attention to the quality of economic growth.Combining with the "12th Five-Year Plan"

Table 2 .
Health index of urban ecosystem in Zibo.

Table 4 .
Main obstacle factors of urban ecosystem health in Zibo, unit: %.